Quantum criticality in the SO(5) bilinear-biquadratic Heisenberg chain 
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The zero-temperature properties of the SO(5) bilinear-biquadratic Heisenberg spin chain are investigated by 
means of a low-energy approach and large scale numerical calculations. In sharp contrast to the spin-1 SO(3) 
Heisenberg chain, we show that the SO(5) Heisenberg spin chain is dimerized with a two-fold degenerate ground 
state. On top of this gapful phase, we find the emergence of a non-degenerate gapped phase with hidden (Z2 x 
Z2) 2 symmetry and spin- 3/2 edge states that can be understood from a SO(5) AKLT wave function. We derive 
a low-energy theory describing the quantum critical point which separates these two gapped phases. It is shown 
and confirmed numerically that this quantum critical point belongs to the SO(5)i universality class. 
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One-dimensional (ID) quantum systems display a wealth 
of fascinating behaviors which have attracted much interest 
over the years. A paradigmatic example of an exotic phase 
stabilized by ID quantum fluctuations is the Haldane phase of 
the spin-1 SU(2) Heisenberg chain fljj] . On top of the existence 
of a gap, this phase displays remarkable properties like the 
existence of a hidden Neel antiferromagnetic order J2] or the 
emergence of fractional spin- 1/2 edge states when the chain 
is doped by non-magnetic impurities (3[]. A simple way to 
grasp the main characteristics of the Haldane phase is pro- 
vided by the seminal work of Affleck, Kennedy, Lieb and- 
Tasaki (AKLT) where a fine-tuned biquadratic exchange 
interaction is introduced so that the ground state (GS) is made 
up solely of nearest-neighbor valence bonds. Due to the ab- 
sence of a quantum phase transition upon switching on the 
biquadratic exchange from the Heisenberg model up to the 
AKLT point, all the exotic properties of the Haldane phase 
can be simply deduced from the AKLT wave function ||5|] . 

In an interesting work, Tu et al. J(J have recently studied 
a generalization of the spin-1 bilinear-biquadratic Heisenberg 
chain to higher symmetry group SO(2rt+l), with Hamiltonian: 
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where L a {, (1 < a < b < 2n + 1) are the n(2n + 1) gener- 
ators which transform in the vectorial representation (n x n 
matrices) of the SO(2rt+l) group. They are normalized such 
that the single-site Casimir operator is: a<b {Lf') 2 = 2n. 
For n = 1, this model is nothing but the spin-1 bilinear bi- 
quadratic Heisenberg chain. The phase diagram of the model 
(Q3 for general n has been conjectured in Ref. 6 based on the 
behaviour at some special points. For 8 — tan -1 l/(2n — 1), 
the model is the SU(2n+l) Sutherland model which displays 
a quantum critical behavior with 2n gapless modes [7]. There 
is a second SU(27i+l) symmetric model for 8 = ±ir/2 with 
alternating fundamental and conjugate representations, where 
the GS is dimerized with a spontaneous breaking of transla- 
tion symmetry yjy. For 8aklt = tan -1 l/(2n + 1), the 



GS is a SO(2n+l) matrix product state which is the gener- 
alization of the SO(5) AKLT state introduced by Scalapino 
et al. [9] in the context of a SO(5) two-leg ladder. The re- 
sulting non-degenerate gapped phase is the SO(2n+l) gener- 
alization of the Haldane phase with hidden antiferromagnetic 
ordering due to the breaking of a (Z2 x Z2)" symmetry 
In open boundary geometry, the two edge states belong to the 
spinorial representation of SO(2n+l) with dimension 2", so 
that the GS is 4 n -fold degenerate. Finally, the model with 
8 R = tan" 1 (27i - 3)/(2n - l) 2 is exactly solvable JT^l and 
is expected to have gapless excitations from general grounds. 

In this paper, we investigate part of the zero-temperature 
phase diagram between 8 = and 8ak lt, including the criti- 
cality of the SO(5) bilinear-biquadratic spin chain by means of 
complementary low-energy approach and several numerical 
calculations: exact diagonalization (ED) and density matrix 
renormalization group (DMRG) ll 111 . We give further argu- 
ments supporting the phase diagram for n — 2 conjectured in 
Ref. la. In sharp contrast with the spin- 1 Heisenberg chain, the 
AKLT point for n = 2 is found not to capture the physics of 
the SO(5) Heisenberg chain. A quantum phase transition oc- 
curs at 8 = Or that we fully characterize and show to belong 
to the SO(5)i universality class with central charge c = 5/2. 

Low-energy approach. In principle, a field-theory approach 
capturing the low-energy properties of the SO(2n+l) Heisen- 
berg spin chain can be derived from the solvable point at 
8 = Or. By exploiting the integrability of the model, one 
may determine the nature of the underlying conformal field 
theory (CFT) and investigate small deviations \6 — 8n\ <C 1 
in parallel to the Majorana fermion approach of Tsvelik for 
n = 1 [ l3- Here, we instead directly derive a low-energy de- 
scription for the SO(5) Heisenberg spin chain by considering a 
spin-3/2 fermionic model which is SO(5) symmetric iflUflill : 
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where i denotes the fermionic creation operator with spin 

index a = ±3/2, ±1/2, and n, = E a c iL i c <M is tne den- 
sity operator at site i. In Eq. d2j, the singlet BCS pair- 
ing operator for spin-3/2 fermions is Pq i — c^ 2 4 c_ 3 , 2 . — 

c i/2 i c -i/2 r ^ s snown i n Refill the spin-3/2 model (f2]l en- 
joys a U(l) c hargc x SO(5) sp in continuous symmetry without 
any fine-tuning. When U, V > 0, the lowest states for t = 
at half-filling are the quintet states: 
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A standard strong-coupling approach can then be applied 
when U, V t > and one finds the effective model to be 
a SO(5) Heisenberg spin chain: H = J£,-£ 0<6 L t L i+i 
with J = 2t 2 /3U. Then, a low-energy description can be 
deduced for the latter model given the adiabatic continuity 
between weak and strong coupling regimes for the spin-3/2 
model [15]. As derived in Ref. Il5l the weak-coupling ap- 
proach at half-filling is built from eight right and left moving 
Majorana fermions £^ L (with A = 1, . . . , 8): 
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where the two Majorana fermions L account for the U(l) 
charge symmetry, the five Majorana fermions ^ L (with a = 
1, . . . , 5) generate the SO(5) symmetry, and ^ R L describes an 
internal discrete Z2 symmetry. By integrating out the charge 
degrees of freedom which are fully gapped in the large U,V 
repulsive limit, one then obtains an effective Hamiltonian with 
leading part 
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describing five massive Majorana fermions with mass m s ~ 
t/ > 0. In fact, after the integration of the charge degrees of 
freedom, £^ L is also a massive Majorana fermion, but with a 
higher mass m > m s so that its contribution can be neglected 
in the low-energy limit E <C too. The effective low-energy 
approach of the SO(5) Heisenberg spin chain is thus given by 
model (SJ, a natural extension to SO(5) of Tsvelik's result (i.e. 
three massive Majorana fermions) for the spin-1 Heisenberg 
chain 1 12]. Since a ID theory of massive Majorana accounts 
for the long-distance behavior of a ID quantum Ising model, 
the low-energy properties of the SO(5) Heisenberg spin chain 
can be described in terms of five decoupled off-critical ID 
quantum Ising models in their disordered phases (m s ~ T — 
T c ). 



To fully characterize the nature of the GS of the SO(5) 
Heisenberg chain, it is instructive to investigate the ex- 
pectation value of the SO(5) dimerization operator Di = 
(— Y E a <b Lf^, One can follow a similar strategy to 
derive a low-energy expression for this operator. In the low- 
energy limit E -C m , we find that D is expressed in terms of 
the disorder operators fi a of the underlying ID quantum Ising 
degrees of freedom: D ~ Ilo=i Since the Ising models 
are locked in their disordered phases, we obtain that (D) ^ 0. 
The GS of the SO(5) Heisenberg spin chain is thus two- 
fold degenerate and spontaneously dimerized in sharp contrast 
with the Haldane phase of the spin-1 Heisenberg chain. 

With this low-energy approach at hand, we can investigate 
the effect of the SO(5) biquadratic term in Eq. (Q3 in the vicin- 
ity of 6 = 0. The most relevant term consistent with the 
symmetries of the lattice model is the Majorana fermion mass 
term of Eq. (0). We thus deduce that the low-energy Hamil- 
tonian of Eq. (Q} is still given by (0]) with a mass m s (9), a 
phenomenological parameter within our approach. The only 
possibility is that the mass m s vanishes and changes sign to 
form a new gapful phase with the five effective ID quan- 
tum Ising models entering their ordered phases. It is very 
tempting to identify the quantum critical point at m s = 
with the location of the integrable model at 9 = 8r which 
is expected, on general grounds, to display a quantum crit- 
ical behavior. This will be confirmed precisely numerically 
below. We thus conclude that the universality class of the 
integrable model for 9 = Or is the SO(5)i CFT with cen- 
tral charge c = 5/2. In addition, our low-energy approach 
enables us to extract the leading asymptotics of the SO(5) 
spin-spin correlations which display a power-law decay with 
a universal exponent: (L^Lf) - (-) x x~ 5 ^. The SO(5) 
dimerization operator vanishes at 9 = 9r and has a power- 
law decaying correlation function with the same universal ex- 
ponent. For 9 > 9r, a new gapful phase emerges with no 
spontaneous dimerization since (// a ) = for m s < 0. This 
phase is the analog of the Haldane phase of the spin-1 Heisen- 
berg chain. Each Ising model has a doubly degenerate GS 
which gives a 2 5 degeneracy. However, there is a redundancy 
in the Majorana fermion description since the transformation 
— > —£,fii J ,H a — > /J, a ,cr a — > —a a leaves invariant the 
Hamiltonian and the physical operators. The correct hidden 
symmetry group is thus (Z2 x Z2) 2 in full agreement with the 
description of the AKLT wave function for the SO(5) general- 
ization of the Haldane phase 10,0]. As a consequence of this 
antiferromagnetic ordering, there are non-trivial edge states 
in this phase when a semi-infinite open chain is considered. 
Using the results of Ref. we find the emergence of five lo- 
calized Majorana fermion zero mode states if 1 inside the gap 
(mid-gap states). These five local fermionic states give rise to 
a local spin-3/2 operator S since: 
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describes a local spin-3/2 operator thanks to the anticommu- 
tation relations of the Majorana fermions: {77°, r/ fc } = S ab . In 
the phase with 9 > 9 b, we thus predict the occurence of spin- 
3/2 boundary excitations at the edge of the chain, again in full 
agreement with the SO(5) AKLT wave-function |9]. 



Numerical results. We combine ED and DMRG calcula- 
tions to clarify the phase diagram and the criticality of the 
phase transition. We use a spin-2 formulation of the model 
(Q} using projectors jfj, which leads to 



H = ^cos^[-l-|s i -S J + (S J -S J ) 2 /9 + (S r S 3 ) 3 /18]+sm0[l-5S J -S,-^(S, 
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(Si-Sj) 4 ]. (5) 



First, we compute with ED the first singlet and magnetic ex- 
citations for various chains with periodic boundary conditions 
(PBC). For 8 < 9r, the data (shown in Fig. [TJ are compatible 
with a finite spin gap but a vanishing singlet gap (with mo- 
mentum 7r), which suggests a singlet phase that breaks trans- 
lation symmetry. It is remarkable that the extrapolation (using 
polynomial fit) of the crossing point of the first magnetic and 
non-magnetic excitations gives 9 C = 6.34° which is precisely 
the numerical value of 9r. 
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FIG. 1: ED data showing the gap to the first singlet and magnetic 
excitations (empty and filled symbols respectively) vs 6 for several 
chain lengths L. Inset shows the extrapolation of the crossing point. 



The dimerized nature of the GS of the Heisenberg SO(5) 
chain at 9 = is confirmed by computing the correla- 
tor (J2 a <b ^i b ^i+i) as a function of position i in an open 
chain. The DMRG data (see Fig. [2]) obtained on a L = 256 
chain with open boundary conditions (OBC) in the range 
9 £ [0, 14°] clearly show a staggered pattern for < 6 < Or, 
in contrast with the uniform values found for 9 > 9r, for 
instance at the AKLT point. 

For 9 ~ 9 aklt, the lowest excitation is magnetic and 
has a finite gap. In fact, right at AKLT point, ED proves 
the existence of a finite gap in the thermodynamic limit us- 
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FIG. 2: DMRG data of the dimerization (J2 a< b L t ' L t+i) asafunc- 
tion of position i, induced in a L = 256 chain by the OBC in the 
ranges 6 < 6r (left panel) and 6 > 6r (right panel). At 8 = 0, the 
system is strongly dimerized. 



ing an argument from Knabe 11711 : the thermodynamic gap is 
bounded by the gap on a finite cluster with L + 1 sites (OBC) 
as Aoo > -j-^y(Al — j;). Using L = 8, we already obtain 
that the spin gap is finite and satisfies Aco > 0.185. 

For larger 9, data are also compatible with a SU(5) gapless 
phase but we observe strong finite-size effects due to incom- 
mensurate excitations, similar to the case of the spin-1 chain 
close to SU(3) point I3. 

As argued in the low-energy section, the central charge of 
the critical point between the dimerized and SO(5) Haldane 
phases should be c = 5/2, corresponding to a SO(5)i CFT In 
Fig. [3^, we show in log-linear scale the von Neumann entropy 
S(x) of a block of x sites at 9r as a function of the conformal 
distance d(x) = L/irsm(irx/L), for chains with PBC of sizes 
up to L = 32. Fitting the data to the expected form for PBC 
S(x) = c/3\n(d(x)) + K with K a constant (see Ref. I20I), 
we obtain an estimate of the central charge c = 2.53. We can 
also obtain c using finite-size corrections of the GS energy per 
site: eo(L) = eo — (irvc)/6L 2 . Using ED data (not shown) 
for both the energy and finite-size velocity v(L) = Eo(k = 
2tt/L) — Eo(k = 0), we obtain the same estimate c = 2.53, 
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in excellent agreement with a SO(5)i quantum criticality. 
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FIG. 3: At 8r and as a function of conformal distance d(x) = 
L/-k sm(x7r /L) for systems with PBC: (a) von Neumann entropy 
of a block of size x (log4in scale), (b) Spin ((— ) E So • S x ) and SO(5) 
generators ((— ) x J2 a <b^o b Lf') correlators (log-log scale). Corre- 
lators are translation-invariant due to the PBC and are normalized to 
1 at their onsite values. 



At the critical point Or, the low-energy approach predicts 
the leading behavior of two-point correlations: (£" Lfi_ x ) ~ 
(_y x -s/4 (the 

same dependence is expected for spin corre- 
lations (Si • Si+a;)). These correlators, plotted in Fig. [3p. for 
a L = 32 chain with PBC, exhibit an even/odd effect besides 
the (— ) x factor. Fits of the two sets of points respectively 
lead to power-law exponents of 1.25 and 1.12, equal or close 
to the theoretical prediction 5/4 for both correlators. We ex- 
pect that all correlators decay with the same exponent 5/4 in 
the thermodynamic limit and that the even/odd effect (and the 
corresponding different exponents) is caused by subleading 
correlations, which can be important for the moderate system 
size L — 32 that can be reached in the DMRG simulations. 

Conclusion. Using complementary techniques, we confirm 
the phase diagram that has been proposed phenomenologi- 
cally in Ref. 6. In particular, we derive a low-energy approach 
in terms of five massive Majorana fermions which allows to 
determine that (i) the GS at 9 = is spontaneously dimer- 
ized and (ii) it is separated from the SO(5) generalization of 
the Haldane phase by a quantum critical point described by a 
SO(5)i CFT with central charge c = 5/2. All these predic- 
tions are quantitatively confirmed using extensive numerical 
simulations. In sharp contrast to the n — 1 case, the SO(5) 
AKLT point does not share the same physics with the Heisen- 
berg point 6 = 0. This property should be a general feature of 
the model with n > 1 |@1 ■ In this respect, the n — 1 case is 
very special and not representative of the generic situation. 

To give some perspective, let us mention that the low- 
energy approach presented here is directly relevant to the 
quantum critical behavior of the spin-2 model interpolating 
between two topologically distinct AKLT states !12lH23ll . In 
particular, the quantum phase transition between the SO(5) 



AKLT state and the spin-2 dimerized phase found in Refs. 

can be shown to be described by the effective Hamil- 
tonian (|4]i and thus belongs to the SO(5)i universality class. 
However, our approach does not capture the direct topologi- 
cal quantum phase transition between the spin-2 AKLT state 
with spin-1 edg e states and the SO(5) AKLT state put forward 
in Refs. I2ll423l It will be interesting to understand how it can 
be extended to explain the emergence of SO(5)i criticality re- 
cently observed numerically for this transition 12311 . 
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